Comparative structural and functional analysis of the glycine-rich regions of Class A and B J-domain protein cochaperones of Hsp70.

J-domain proteins are critical Hsp70 co-chaperones. A and B types have a poorly understood glycine-rich region (Grich) adjacent to their N-terminal J-domain (Jdom). We analyzed the ability of Jdom/Grich segments of yeast Class B Sis1 and a suppressor variant of Class A, Ydj1, to rescue the inviability of sis1-∆. In each, we identified a cluster of Grich residues required for rescue. Both contain conserved hydrophobic and acidic residues and are predicted to form helices. While, as expected, the Sis1 segment docks on its J-domain, that of Ydj1 does not. However, data suggest both interact with Hsp70. We speculate that the Grich-Hsp70 interaction of Classes A and B J-domain proteins can fine tune the activity of Hsp70, thus being particularly important for the function of Class B.

J-domain proteins: From molecular mechanisms to diseases.

J-domain proteins (JDPs) are the largest family of chaperones in most organisms, but much of how they function within the network of other chaperones and protein quality control machineries is still an enigma. Here, we report on the latest findings related to JDP functions presented at a dedicated JDP workshop in Gdansk, Poland. The report does not include all (details) of what was shared and discussed at the meeting, because some of these original data have not yet been accepted for publication elsewhere or represented still preliminary observations at the time.

Analysis of Reconstituted Tripartite Complex Supports Avidity-based Recruitment of Hsp70 by Substrate Bound J-domain Protein.

Hsp70 are ubiquitous, versatile molecular chaperones that cyclically interact with substrate protein(s). The initial step requires synergistic interaction of a substrate and a J-domain protein (JDP) cochaperone, via its J-domain, with Hsp70 to stimulate hydrolysis of its bound ATP. This hydrolysis drives conformational changes in Hsp70 that stabilize substrate binding. However, because of the transient nature of substrate and JDP interactions, this key step is not well understood. Here we leverage a well characterized Hsp70 system specialized for iron-sulfur cluster biogenesis, which like many systems, has a JDP that binds substrate on its own. Utilizing an ATPase-deficient Hsp70 variant, we isolated a Hsp70-JDP-substrate tripartite complex. Complex formation and stability depended on residues previously identified as essential for bipartite interactions: JDP-substrate, Hsp70-substrate and J-domain-Hsp70. Computational docking based on the established J-domain-Hsp70(ATP) interaction placed the substrate close to its predicted position in the peptide-binding cleft, with the JDP having the same architecture as when in a bipartite complex with substrate. Together, our results indicate that the structurally rigid JDP-substrate complex recruits Hsp70(ATP) via precise positioning of J-domain and substrate at their respective interaction sites - resulting in functionally high affinity (i.e., avidity). The exceptionally high avidity observed for this specialized system may be unusual because of the rigid architecture of its JDP and the additional JDP-Hsp70 interaction site uncovered in this study. However, functionally important avidity driven by JDP-substrate interactions is likely sufficient to explain synergistic ATPase stimulation and efficient substrate trapping in many Hsp70 systems.

J-Domain Proteins Orchestrate the Multifunctionality of Hsp70s in Mitochondria: Insights from Mechanistic and Evolutionary Analyses.

Mitochondrial J-domain protein (JDP) co-chaperones orchestrate the function of their Hsp70 chaperone partner(s) in critical organellar processes that are essential for cell function. These include folding, refolding, and import of mitochondrial proteins, maintenance of mitochondrial DNA, and biogenesis of iron-sulfur cluster(s) (FeS), prosthetic groups needed for function of mitochondrial and cytosolic proteins. Consistent with the organelle's endosymbiotic origin, mitochondrial Hsp70 and the JDPs' functioning in protein folding and FeS biogenesis clearly descended from bacteria, while the origin of the JDP involved in protein import is less evident. Regardless of their origin, all mitochondrial JDP/Hsp70 systems evolved unique features that allowed them to perform mitochondria-specific functions. Their modes of functional diversification and specialization illustrate the versatility of JDP/Hsp70 systems and inform our understanding of system functioning in other cellular compartments.

Two Bacterial Small Heat Shock Proteins, IbpA and IbpB, Form a Functional Heterodimer.

Small heat shock proteins (sHsps) are a conserved class of ATP-independent chaperones which in stress conditions bind to unfolded protein substrates and prevent their irreversible aggregation. Substrates trapped in sHsps-containing aggregates are efficiently refolded into native structures by ATP-dependent Hsp70 and Hsp100 chaperones. Most γ-proteobacteria possess a single sHsp (IbpA), while in a subset of Enterobacterales, as a consequence of ibpA gene duplication event, a two-protein sHsp (IbpA and IbpB) system has evolved. IbpA and IbpB are functionally divergent. Purified IbpA, but not IbpB, stably interacts with aggregated substrates, yet both sHsps are required to be present at the substrate denaturation step for subsequent efficient Hsp70-Hsp100-dependent substrate refolding. IbpA and IbpB interact with each other, influence each other's expression levels and degradation rates. However, the crucial information on how these two sHsps interact and what is the basic building block required for proper sHsps functioning was missing. Here, based on NMR, mass spectrometry and crosslinking studies, we show that IbpA-IbpB heterodimer is a dominating functional unit of the two sHsp system in Enterobacterales. The principle of heterodimer formation is similar to one described for homodimers of single bacterial sHsps. ß-hairpins formed by strands ß5 and ß7 of IbpA or IbpB crystallin domains associate with the other one's ß-sandwich in the heterodimer structure. Relying on crosslinking and molecular dynamics studies, we also propose the orientation of two IbpA-IbpB heterodimers in a higher order tetrameric structure.

Defining a novel domain that provides an essential contribution to site-specific interaction of Rep protein with DNA.

An essential feature of replication initiation proteins is their ability to bind to DNA. In this work, we describe a new domain that contributes to a replication initiator sequence-specific interaction with DNA. Applying biochemical assays and structure prediction methods coupled with DNA-protein crosslinking, mass spectrometry, and construction and analysis of mutant proteins, we identified that the replication initiator of the broad host range plasmid RK2, in addition to two winged helix domains, contains a third DNA-binding domain. The phylogenetic analysis revealed that the composition of this unique domain is typical within the described TrfA-like protein family. Both in vitro and in vivo experiments involving the constructed TrfA mutant proteins showed that the newly identified domain is essential for the formation of the protein complex with DNA, contributes to the avidity for interaction with DNA, and the replication activity of the initiator. The analysis of mutant proteins, each containing a single substitution, showed that each of the three domains composing TrfA is essential for the formation of the protein complex with DNA. Furthermore, the new domain, along with the winged helix domains, contributes to the sequence specificity of replication initiator interaction within the plasmid replication origin.

Two-step mechanism of J-domain action in driving Hsp70 function.

J-domain proteins (JDPs), obligatory Hsp70 cochaperones, play critical roles in protein homeostasis. They promote key allosteric transitions that stabilize Hsp70 interaction with substrate polypeptides upon hydrolysis of its bound ATP. Although a recent crystal structure revealed the physical mode of interaction between a J-domain and an Hsp70, the structural and dynamic consequences of J-domain action once bound and how Hsp70s discriminate among its multiple JDP partners remain enigmatic. We combined free energy simulations, biochemical assays and evolutionary analyses to address these issues. Our results indicate that the invariant aspartate of the J-domain perturbs a conserved intramolecular Hsp70 network of contacts that crosses domains. This perturbation leads to destabilization of the domain-domain interface-thereby promoting the allosteric transition that triggers ATP hydrolysis. While this mechanistic step is driven by conserved residues, evolutionarily variable residues are key to initial JDP/Hsp70 recognition-via electrostatic interactions between oppositely charged surfaces. We speculate that these variable residues allow an Hsp70 to discriminate amongst JDP partners, as many of them have coevolved. Together, our data points to a two-step mode of J-domain action, a recognition stage followed by a mechanistic stage.

Duplicate divergence of two bacterial small heat shock proteins reduces the demand for Hsp70 in refolding of substrates.

Small heat shock proteins (sHsps) are a conserved class of ATP-independent chaperones that bind to aggregation-prone polypeptides at stress conditions. sHsps encage these polypeptides in assemblies, shielding them from further aggregation. To facilitate their subsequent solubilization and refolding by Hsp70 (DnaK) and Hsp100 (ClpB) chaperones, first, sHsps need to dissociate from the assemblies. In most γ-proteobacteria, these functions are fulfilled by a single sHsp (IbpA), but in a subset of Enterobacterales, a two-protein sHsp (IbpA and IbpB) system has evolved. To gain insight into the emergence of complexity within this chaperone system, we reconstructed the phylogeny of γ-proteobacteria and their sHsps. We selected proteins representative of systems comprising either one or two sHsps and analysed their ability to form sHsps-substrate assemblies. All the tested IbpA proteins, but not IbpBs, stably interact with an aggregating substrate. Moreover, in Escherichia coli cells, ibpA but not ibpB suppress the growth defect associated with low DnaK level, which points to the major protective role of IbpA during the breakdown of protein quality control. We also examined how sHsps affect the association of Hsp70 with the assemblies at the initial phase of disaggregation and how they affect protein recovery after stress. Our results suggest that a single gene duplication event has given rise to the sHsp system consisting of a strong canonical binder, IbpA, and its non-canonical paralog IbpB that enhances sHsps dissociation from the assemblies. The cooperation between the sHsps reduces the demand for Hsp70 needed to outcompete them from the assemblies by promoting sHsps dissociation without compromising assembly formation at heat shock. This potentially increases the robustness and elasticity of sHsps protection against irreversible aggregation.

OMA standalone: orthology inference among public and custom genomes and transcriptomes.

Genomes and transcriptomes are now typically sequenced by individual laboratories but analyzing them often remains challenging. One essential step in many analyses lies in identifying orthologs-corresponding genes across multiple species-but this is far from trivial. The Orthologous MAtrix (OMA) database is a leading resource for identifying orthologs among publicly available, complete genomes. Here, we describe the OMA pipeline available as a standalone program for Linux and Mac. When run on a cluster, it has native support for the LSF, SGE, PBS Pro, and Slurm job schedulers and can scale up to thousands of parallel processes. Another key feature of OMA standalone is that users can combine their own data with existing public data by exporting genomes and precomputed alignments from the OMA database, which currently contains over 2100 complete genomes. We compare OMA standalone to other methods in the context of phylogenetic tree inference, by inferring a phylogeny of Lophotrochozoa, a challenging clade within the protostomes. We also discuss other potential applications of OMA standalone, including identifying gene families having undergone duplications/losses in specific clades, and identifying potential drug targets in nonmodel organisms. OMA standalone is available under the permissive open source Mozilla Public License Version 2.0.

Structure and evolution of the 4-helix bundle domain of Zuotin, a J-domain protein co-chaperone of Hsp70.

The J-domain protein Zuotin is a multi-domain eukaryotic Hsp70 co-chaperone. Though it is primarily ribosome-associated, positioned at the exit of the 60S subunit tunnel where it promotes folding of nascent polypeptide chains, Zuotin also has off-ribosome functions. Domains of Zuotin needed for 60S association and interaction with Hsp70 are conserved in eukaryotes. However, whether the 4-helix bundle (4HB) domain is conserved remains an open question. We undertook evolutionary and structural approaches to clarify this issue. We found that the 4HB segment of human Zuotin also forms a bundle of 4 helices. The positive charge of Helix I, which in Saccharomyces cerevisiae is responsible for interaction with the 40S subunit, is particularly conserved. However, the C-termini of fungal and human 4HBs are not similar. In fungi the C-terminal segment forms a plug that folds back into the bundle; in S. cerevisiae it plays an important role in bundle stability and, off the ribosome, in transcriptional activation. In human, C-terminal helix IV of the 4HB is extended, protruding from the bundle. This extension serves as a linker to the regulatory SANT domains, which are present in animals, plants and protists, but not fungi. Further analysis of Zuotin sequences revealed that the plug likely arose as a result of genomic rearrangement upon SANT domain loss early in the fungal lineage. In the lineage leading to S. cerevisiae, the 4HB was subjected to positive selection with the plug becoming increasingly hydrophobic. Eventually, these hydrophobic plug residues were coopted for a novel regulatory function-activation of a recently emerged transcription factor, Pdr1. Our data suggests that Zuotin evolved off-ribosome functions twice-once involving SANT domains, then later in fungi, after SANT domain loss, by coopting the hydrophobic plug. Zuotin serves as an example of complex intertwining of molecular chaperone function and cell regulation.

Mitigating Anticipated Effects of Systematic Errors Supports Sister-Group Relationship between Xenacoelomorpha and Ambulacraria.

Xenoturbella and the acoelomorph worms (Xenacoelomorpha) are simple marine animals with controversial affinities. They have been placed as the sister group of all other bilaterian animals (Nephrozoa hypothesis), implying their simplicity is an ancient characteristic [1, 2]; alternatively, they have been linked to the complex Ambulacraria (echinoderms and hemichordates) in a clade called the Xenambulacraria [3-5], suggesting their simplicity evolved by reduction from a complex ancestor. The difficulty resolving this problem implies the phylogenetic signal supporting the correct solution is weak and affected by inadequate modeling, creating a misleading non-phylogenetic signal. The idea that the Nephrozoa hypothesis might be an artifact is prompted by the faster molecular evolutionary rate observed within the Acoelomorpha. Unequal rates of evolution are known to result in the systematic artifact of long branch attraction, which would be predicted to result in an attraction between long-branch acoelomorphs and the outgroup, pulling them toward the root [6]. Other biases inadequately accommodated by the models used can also have strong effects, exacerbated in the context of short internal branches and long terminal branches [7]. We have assembled a large and informative dataset to address this problem. Analyses designed to reduce or to emphasize misleading signals show the Nephrozoa hypothesis is supported under conditions expected to exacerbate errors, and the Xenambulacraria hypothesis is preferred in conditions designed to reduce errors. Our reanalyses of two other recently published datasets [1, 2] produce the same result. We conclude that the Xenacoelomorpha are simplified relatives of the Ambulacraria.

Iron-Sulfur Cluster Biogenesis Chaperones: Evidence for Emergence of Mutational Robustness of a Highly Specific Protein-Protein Interaction.

Biogenesis of iron-sulfur clusters (FeS) is a highly conserved process involving Hsp70 and J-protein chaperones. However, Hsp70 specialization differs among species. In most eukaryotes, including Schizosaccharomyces pombe, FeS biogenesis involves interaction between the J-protein Jac1 and the multifunctional Hsp70 Ssc1. But, in Saccharomyces cerevisiae and closely related species, Jac1 interacts with the specialized Hsp70 Ssq1, which emerged through duplication of SSC1. As little is known about how gene duplicates affect the robustness of their protein interaction partners, we analyzed the functional and evolutionary consequences of Ssq1 specialization on the ubiquitous J-protein cochaperone Jac1, by comparing S. cerevisiae and S. pombe. Although deletion of JAC1 is lethal in both species, alanine substitutions within the conserved His-Pro-Asp (HPD) motif, which is critical for Jac1:Hsp70 interaction, have species-specific effects. They are lethal in S. pombe, but not in S. cerevisiae. These in vivo differences correlated with in vitro biochemical measurements. Charged residues present in the J-domain of S. cerevisiae Jac1, but absent in S. pombe Jac1, are important for tolerance of S. cerevisiae Jac1 to HPD alterations. Moreover, Jac1 orthologs from species that encode Ssq1 have a higher sequence divergence. The simplest interpretation of our results is that Ssq1's coevolution with Jac1 resulted in expansion of their binding interface, thus increasing the efficiency of their interaction. Such an expansion could in turn compensate for negative effects of HPD substitutions. Thus, our results support the idea that the robustness of Jac1 emerged as consequence of its highly efficient and specific interaction with Ssq1.

A transcriptomic-phylogenomic analysis of the evolutionary relationships of flatworms.

The interrelationships of the flatworms (phylum Platyhelminthes) are poorly resolved despite decades of morphological and molecular phylogenetic studies. The earliest-branching clades (Catenulida, Macrostomorpha, and Polycladida) share spiral cleavage and entolecithal eggs with other lophotrochozoans. Lecithoepitheliata have primitive spiral cleavage but derived ectolecithal eggs. Other orders (Rhabdocoela, Proseriata, Tricladida and relatives, and Bothrioplanida) all have derived ectolecithal eggs but have uncertain affinities to one another. The orders of parasitic Neodermata emerge from an uncertain position from within these ectolecithal classes. To tackle these problems, we have sequenced transcriptomes from 18 flatworms and 5 other metazoan groups. The addition of published data produces an alignment of >107,000 amino acids with less than 28% missing data from 27 flatworm taxa in 11 orders covering all major clades. Our phylogenetic analyses show that Platyhelminthes consist of the two clades Catenulida and Rhabditophora. Within Rhabditophora, we show the earliest-emerging branch is Macrostomorpha, not Polycladida. We show Lecithoepitheliata are not members of Neoophora but are sister group of Polycladida, implying independent origins of the ectolecithal eggs found in Lecithoepitheliata and Neoophora. We resolve Rhabdocoela as the most basally branching euneoophoran taxon. Tricladida, Bothrioplanida, and Neodermata constitute a group that appears to have lost both spiral cleavage and centrosomes. We identify Bothrioplanida as the long-sought closest free-living sister group of the parasitic Neodermata. Among parasitic orders, we show that Cestoda are closer to Trematoda than to Monogenea, rejecting the concept of the Cercomeromorpha. Our results have important implications for understanding the evolution of this major phylum.

The genomic and phenotypic diversity of Schizosaccharomyces pombe.

Natural variation within species reveals aspects of genome evolution and function. The fission yeast Schizosaccharomyces pombe is an important model for eukaryotic biology, but researchers typically use one standard laboratory strain. To extend the usefulness of this model, we surveyed the genomic and phenotypic variation in 161 natural isolates. We sequenced the genomes of all strains, finding moderate genetic diversity (π = 3 × 10(-3) substitutions/site) and weak global population structure. We estimate that dispersal of S. pombe began during human antiquity (∼340 BCE), and ancestors of these strains reached the Americas at ∼1623 CE. We quantified 74 traits, finding substantial heritable phenotypic diversity. We conducted 223 genome-wide association studies, with 89 traits showing at least one association. The most significant variant for each trait explained 22% of the phenotypic variance on average, with indels having larger effects than SNPs. This analysis represents a rich resource to examine genotype-phenotype relationships in a tractable model.

The OMA orthology database in 2015: function predictions, better plant support, synteny view and other improvements.

The Orthologous Matrix (OMA) project is a method and associated database inferring evolutionary relationships amongst currently 1706 complete proteomes (i.e. the protein sequence associated for every protein-coding gene in all genomes). In this update article, we present six major new developments in OMA: (i) a new web interface; (ii) Gene Ontology function predictions as part of the OMA pipeline; (iii) better support for plant genomes and in particular homeologs in the wheat genome; (iv) a new synteny viewer providing the genomic context of orthologs; (v) statically computed hierarchical orthologous groups subsets downloadable in OrthoXML format; and (vi) possibility to export parts of the all-against-all computations and to combine them with custom data for 'client-side' orthology prediction. OMA can be accessed through the OMA Browser and various programmatic interfaces at http://omabrowser.org.

A beginners guide to estimating the non-synonymous to synonymous rate ratio of all protein-coding genes in a genome.

The ratio of non-synonymous to synonymous substitutions (dN/dS) is a useful measure of the strength and mode of natural selection acting on protein-coding genes. It is widely used to study patterns of selection on protein genes on a genomic scale-from the small genomes of viruses, bacteria, and parasitic eukaryotes to the largest eukaryotic genomes. In this chapter we describe all the steps necessary to calculate the dN/dS of all the genes using at least two genomes. We include a brief discussion on assigning orthologs, and of codon-aware alignment of orthologs. We then describe how to use the CODEML program of the PAML package for phylogenetic analysis to calculate the dN/dS and how to perform some statistical tests for positive selection. We then outline some methods for interpreting output and describe how one may use this data to make discoveries about the biology of your species. Finally, as a worked example we show all the steps we used to calculate dN/dS for 3,261 orthologs from six Plasmodium species, including tests for adaptive evolution (see worked_example.pdf).

Molecular recognition in complexes of TRF proteins with telomeric DNA.

Telomeres are specialized nucleoprotein assemblies that protect the ends of linear chromosomes. In humans and many other species, telomeres consist of tandem TTAGGG repeats bound by a protein complex known as shelterin that remodels telomeric DNA into a protective loop structure and regulates telomere homeostasis. Shelterin recognizes telomeric repeats through its two major components known as Telomere Repeat-Binding Factors, TRF1 and TRF2. These two homologous proteins are therefore essential for the formation and normal function of telomeres. Indeed, TRF1 and TRF2 are implicated in a plethora of different cellular functions and their depletion leads to telomere dysfunction with chromosomal fusions, followed by apoptotic cell death. More specifically, it was found that TRF1 acts as a negative regulator of telomere length, and TRF2 is involved in stabilizing the loop structure. Consequently, these proteins are of great interest, not only because of their key role in telomere maintenance and stability, but also as potential drug targets. In the current study, we investigated the molecular basis of telomeric sequence recognition by TRF1 and TRF2 and their DNA binding mechanism. We used molecular dynamics (MD) to calculate the free energy profiles for binding of TRFs to telomeric DNA. We found that the predicted binding free energies were in good agreement with experimental data. Further, different molecular determinants of binding, such as binding enthalpies and entropies, the hydrogen bonding pattern and changes in surface area, were analyzed to decompose and examine the overall binding free energies at the structural level. With this approach, we were able to draw conclusions regarding the consecutive stages of sequence-specific association, and propose a novel aspartate-dependent mechanism of sequence recognition. Finally, our work demonstrates the applicability of computational MD-based methods to studying protein-DNA interactions.

Put a tiger in your tank: the polyclad flatworm Maritigrella crozieri as a proposed model for evo-devo.

Polyclad flatworms are an early branching clade within the rhabditophoran Platyhelminthes. They provide an interesting system with which to explore the evolution of development within Platyhelminthes and amongst Spiralia (Lophotrochozoa). Unlike most other flatworms, polyclads undergo spiral cleavage (similar to that seen in some other spiralian taxa), they are the only free-living flatworms where development via a larval stage occurs, and they are the only flatworms in which embryos can be reared outside of their protective egg case, enabling embryonic manipulations. Past work has focused on comparing early cleavage patterns and larval anatomy between polyclads and other spiralians. We have selected Maritigrella crozieri, the tiger flatworm, as a suitable polyclad species for developmental studies, because it is abundant and large in size compared to other species. These characteristics have facilitated the generation of a transcriptome from embryonic and larval material and are enabling us to develop methods for gene expression analysis and immunofluorescence techniques. Here we give an overview of M. crozieri and its development, we highlight the advantages and current limitations of this animal as a potential evo-devo model and discuss current lines of research.

Proteasome inhibition can impair caspase-8 activation upon submaximal stimulation of apoptotic tumor necrosis factor-related apoptosis inducing ligand (TRAIL) signaling.

Tumor necrosis factor-related apoptosis inducing ligand (TRAIL) can induce extrinsic apoptosis, resulting in caspase-8 activation, but may also initiate transcription-dependent prosurvival signaling. Proteasome inhibitors were suggested to promote TRAIL signal transduction through the death-inducing signaling complex (DISC) by modulating the relative abundance of core DISC components, thereby enhancing caspase-8 activation and apoptosis. To test this hypothesis, we quantified the changes in DISC protein levels as an early consequence of proteasome inhibition in HeLa cervical cancer cells and, based on these data, mathematically modeled the proapoptotic TRAIL signaling toward caspase-8 activation. Modeling results surprisingly suggested that caspase-8 activation might be delayed in presence of proteasome inhibitors, in particular at submaximal TRAIL doses. Subsequent FRET-based single cell time-lapse imaging at conditions where transcription dependent prosurvival signaling was blocked confirmed this hypothesis: caspase-8 activity was delayed by hours in the presence of proteasome inhibitors epoxomicin or bortezomib. Corresponding delays were detected for effector caspase processing and cell death. Contrary to current models, we therefore provide evidence that synergies between TRAIL and proteasome inhibitors do not result from changes in the levels of core DISC signaling proteins.